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ABSTRACT 
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Context. M87 is a giant elliptical galaxy located in the centre of the Virgo cluster, which harbours a supermassive black hole of 
mass 6.4 x 10^ whose activity is responsible for the extended (80 kpc) radio lobes that surround the galaxy. The energy 

generated by matter falling onto the central black hole is ejected and transferred to the intra-cluster medium via a relativist ic jet 
and morphologically complex systems of buoyant bubbles, which rise towards the edges of the extended halo. 

Aims. In order to place constraints on past activity cycles of the active nucleus, images of M87 were produced at low radio 
frequencies never explored before at these high spatial resolution and dynamic range. To disentangle different synchrotron models 
and place constraints on source magnetic field, age and energetics, we also performed a detailed spectral analysis of M87 extended 
i radio- halo. 

' Methods. Here we present the first observations made with the new Low-Frequency Array (LOFAR) of M87 at frequencies down 
^ ' to 20 MHz. Three observations were conducted, at 15 — 30 MHz, 30 — 77 MHz and 116 — 162 MHz. We used these observations 
together with archival data to produce a low-frequency spectral index map and to perform a spectral analysis in the wide frequency 
range 30 MHz - 10 GHz. 

■ Results. We do not find any sign of new extended emissions; on the contrary the source appears well confined by the high pressure 
of the intra-cluster medium. A continuous injection of relativistic electrons is the model that best fits our data, and provides a 
scenario in which the lobes are still supplied by fresh relativistic particles from the active galactic nuclei. We suggest that the 
discrepancy between the low-frequency radio-spectral slope in the core and in the halo implies a strong adiabatic expansion of 
the plasma as soon as it leaves the core area. The extended halo has an equipartition magnetic field strength of :^ 10 /iG, which 
increases to 13 /iG in the zones where the particle flows are more active. The continuous injection model for synchrotron ageing 
provides an age for the halo of ^ 40 Myr, which in turn provides a jet kinetic power of 6 — 10 x 10^^ erg s~^. 

Key words, radiation mechanisms: non-thermal - galaxies: active - galaxies: clusters: individual (Virgo) - galaxies: individual (M87) 
- galaxies: jets - radio continuum: galaxies 

1. Introduction amples of black hole - host galaxy feedback in action is the 

AGN in the nearby giant elliptical galaxy M87 (NGC 4486), 

Accreting supermassive black holes in active galactic nuclei -^^ ^^^^ Virgo cluster 

(AGN) can release enormous amounts of energy into their . . i i i r i • i t i 

surroundings, which may profoundly influence the black , ^ particularly large amount of study, mcludmg hun- 

hole's hosting environment up to the cluster scale. The ^^^^^ of published papers, has been devoted to M87. This 

way energy is transported to large distances, its amount §^1^^^ ^^^.^ its popularjrty to several reasons, among others: 

and the typical time-scales of these processes are still not ^V', ^^^^fU radio galaxies it is at the centre 

clear. New-generation radio-telescopes such as the Low- of the nearest rich cluster of galaxies (the Virgo cluster), 

FVequency Array (LOFAR) allow us to study these mecha- ^^^f fo^'^th brightest radio source m the northern sky, 

nisms with unprecedented quality and resolution in a hith- it hosts m its nucleus one of the most massive active 

erto neglected wavelength range. One of the best studied ex- 
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M87 is at a distance of ^ 17 Mpc, where l" corresponds to 
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black holes discovered so f ar (Mbh — 6.4 ± 0.5 x 10^ Mq, 
iGebhardt Thomasl[2QQ9l) . 

The interaction between the AGN of M87 with its 
host galaxy and the intra-cluster medium (ICM) has 
been the subject of a large fraction of the aforemen- 
tioned studies. The emission generated directly and in- 
direc tly by the AGN has been wi d ely observed at ra- 
dio ( i Boltoi et all [1949"; "Mmsl ll952yBaade & Minkowski 
119541: lOwen et all [2000). infrared (IShi et al. 2007), opti- 
cal (Biretta et al. 1999) and X-r ay (Fabricant et al. 1980i 
Feigelson et al. 1987: Bohringer et al. 1995: Yo ung et all 
20021 : iForman et al.ll2007l : iMillion et^ r2010) wavelengths. 



Theoretical and numerical models to inter pret these obser- 
vatio ns ha ve also be en developed by e.g. IChurazov et al.l 
([200 If ) and lBriiggeiTet al.l (|2002i) . 

The radio source associated with this galaxy is named 
Virgo A (3C 274). Its inner region (1.3' x 0.5') contains a 
collimated relativistic jet, which points towards the north- 
west and is embedded in a halo with a diameter up to 15' 

8 kpc) . Th e extended radio eni i ssion, discovered by 
19521 ) and'Baade & Minkowski (1954), is responsible 
for much of the radio flux, especially at the lower frequen- 
cies. Due to the high surface brightness of the compact 
central region and the relatively faint surface brightness 
of the extended emission, high dynamic-range imaging of 
Vir^o A has always been a big challenge. In the past years 
I Owen et al.l (|2Q 00) presented a high-resolution (7"), high- 
dynamic range map of the halo of M87 observed at 327 MHz 
with the Very Large A rray (VLA). At higher frequencies 
iRottmann et al.l (| 19961 ) mapped the extended Virgo halo 
at 10.55 GHz with the single-dish Effelsberg radio telescope 
at 69" resolution. At lower frequency (74 MHz) a 20" res - 
olution map of Virgo A was made by iKassim et all (I1993D . 
This paper will extend the high resolution imaging into the 
previously almost unexplored very low frequency range of 
15 — 162 MHz and present some of the highest-dynamic- 
range images ever made at these frequencies of extended 
source structures. 

M87 lies at the centre of the Virgo cluster X-ray 
luminous atmo sphere, first de t ected with the Einstein 
Observatory by iFabricant et al.l (|l980l ) . An asymmetry in 
the X-ray emission, in the form of two spectacular outflow- 
like structures extending from the nucleus towards the eas t 
and south-west, was discovered bv iFeigelson et all (|1987[ ). 
who also found a correlation between X-ray and radio emit- 
ting features. One of the first explanations for such a cor- 
relation was that the relativistic electrons that produce 
the synchrotron radio emission were also responsible for 
the inverse Compton scattering of cosmic microwave back- 
ground (CMB) photo ns, thus produci ng X-ray radiation 
(IFeigelson et al.l I1987D . However, Bohri nger et al.l (I1995D 
showed with a ROSAT PSPC observation, that the excess 
emission had a thermal spectrum and it is colder than the 
ambient gas which is at a temperature o f > 2 keV. This 



feature was explained by 'Churazov et al.l (|200lf ) as buoy- 
ant bubbles of cosmic rays, injected into the inner halo (or 
"cocoon") by the relativistic jet, which subsequently rise 
through the cooHng gas at about half the sound speed. 
During their rise they uplift gas at the temperature of 
~ 1 keV from the central regions. Sub-arcsecond Chandra 
X-ray images (Million et al. 2010) confirmed this picture 
and provided an unprecedented view of the physical and 
chemical properties of the ICM. 



Although Virgo A is a unique object because of its prop- 
erties, close proximity and sheer quantity of available data, 
it remains a fundamental example to study the more general 
behaviour of AGNs located at the centre of galaxy clusters. 
Knowledge of its energetics and of the interaction between 
its jets and the ICM, may help solve open problems such 
as the suppr ession of the cooHng flows (for a review see 
iPeterson fc F abian 2006) and the AGN duty cycle. This, in 
turn, will provide important clues on the physical nature 
of AGN feedback in massive g alaxies and on its r e levance 
to a cosmological framework (|Croton et al.l 120061 iFabianI 
[2012). Furthermore, it has been claimed that jet sources 
like Virgo A and its southern sibling Centaurus A are po- 
tential candidate sources for t he production of ultra-high 
energ y cosmic rays (UHECR, H ierre Auger Collaboration! 
I2007D . 

In this paper we extend the study of Virgo A to long, 
so far unobserved, wavelengths. We also retrieved available 
observations of Virgo A at 1.4 and 1.6 GHz (VLA, from the 
data archive), at 325 MHz (VLA, provided by Frazer Owen) 
and at 10.55 GHz (Effelsberg radio telescope, provided by 
Helge Rottmann). This enabled us to assess the source en- 
ergetics, the halo age and the main mechanisms which con- 
tributed to its spectral evolution. The paper is organized as 
follows: in the next section we outline the LOFAR features 
and characteristics. In Sect. [3] we present new LOFAR ob- 
servations of Virgo A and we describe the data reduction 
technique. In Sects. H] and [5] we respectively present the 
outcome of these observations and perform a spectral anal- 
ysis of them, discussing the physical interpretation of our 
results. In Sects. [6] and [71 we discuss the results and outline 
our conclusions. 



2. LOFAR 

LOFAR (van Haarlem et al. in prep.) is a radio telescope 
optimized for the frequency range from 30 to 240 MHz, but 
also with the ability to observe down to 10 MHz. LOFAR 
does not have any moving parts, the telescope receivers 
are two different kind of dipoles: the low-band antennas 
(LB A), which cover the frequency range 10 — 90 MHz, and 
the high-band antennas (HBA), which cover the frequency 
range 110 — 240 MHz. The LB As are inverted- V crossed- 
dipoles oriented NE-SW and SE-NW, while the HBA are 
organized into tiles made of a 4 x 4 array of bowtie-shaped 
crossed dipoles. Dipoles are organized into stations yield- 
ing, for each station, effective aperture sizes that range from 
30 m to 80 m, depending on the frequency. Each set of 
dipoles within a station works as a phased aperture array 
- i.e. a delay is applied to the relative phases of the sig- 
nals feeding the dipoles in such a way that the radiation 
pattern of the array is reinforced in a target direction and 
suppressed in undesired ones. By applying different delays, 
LOFAR can therefore "point" (create a beam) in more than 
one direction simultaneously and the number of beams is 
limited only by the bandwidth necessary to transfer the 
signal to the correlator and its computational power. 

The complete configuratior|3 of LOFAR will consist of 
an array of stations distributed over 100 km within the 
Netherlands and out to 1000 km throughout Europe, which 
will provide sufficient resolution \" at 30 MHz) to allow 



^ An updated map of the station status can be f ound here: 
jhttp://www.astron.nl/^heald/lofarSt atusMap.html 
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Fig. 1: 'i/i;-coverage for the three observations of Virgo A: the first row is the HBA observation, the second row the 
LBA-high observation and the third row the LBA-low observation. In the first column are plotted only tracks involving 
remote stations (blue: core-remote baselines - red: remote-remote baselines). In the second column only tracks of core-core 
baselines are plotted. The last column is a zoom-in on the centre of the lit'-plane. 



optical identification of radio sources, even at low frequen- 
cies. At present, international stations provide LOFAR with 
an angular resolution of ~ 0".15 at 240 MHz and ^ at 
30 MHz, while a dense core of 24 stations provides the nec- 
essary sensitivity to the extended emission. 



For the stations in the Netherlands only 48 (out of 
96) LBA dipoles can be currently used simultaneously. 
Among the different possibilities, mainly two possible con- 
figurations of the LBA dipoles are commonly used: LBA- 
INNER and LBA-OUTER. In these configurations the ac- 
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Fig. 2: Left: dirty beam used in the deconvolution step for the HBA observation (at a frequency of ~ 140 MHz). Centre: 
dirty beam used in the deconvolution step for the LB A- high observation (at a frequency of ~ 54 MHz). Right: dirty 
beam used in the deconvolution step for the LB A- low observation (at a frequency of ^ 22 MHz). 



tive dipoles are located respectively in the inner zone and in 
the outer zone of the station field. The more concentrated 
are the dipoles in the inner zone, the more the side-lobe 
levels are reduced and the field of view (FoV) is wider, but 
at the cost of a reduced sensitivity. Each core station has 
the HBA grouped into two different sub-stations located at 
the edge of the field. These sub-stations can be used to- 
gether or as standalone stations (DUAL observing mode), 
to increase the number of baselines. 

Once the data are collected from the stations they 
are transported to the central processing location via a 
Wide- Area Network, using dedicated light paths. Data are 
then correlated by a Blue Gene/P computer that contains 
12480 processor cores providing 42.4 TFLOPS peak pro- 
ces sing power. Fo r a de tailed description of the correlator 
see lRomein et al.l (|2010[ ). 



3. The observations 

In this paper we present a set of three observations per- 
formed during the LOFAR commissioning phase. The phase 
centre was set on the core of Virgo A (RA: 12:30:49.420 - 
DEC: +12:23:28.0 - J2000) and the observational details 
are listed in Table [H Each observation was 8 hours in dura- 
tion and all four polarization products (XX, YY, XY, and 
YX) were stored. Each observation had its frequency cov- 
erage divided into sub-bands (SB) of 0.1953 MHz of band- 
width and each SB was divided into 64 channels of 3 kHz 
of bandwidth. The following observations were performed: 

HBA (115 - 162 MHz): we observed the target with the 
HBA on the 2nd and 3rd of April, 2011. The visibil- 
ity sampling time was 2 s. Two stations (CS021HBA0 
and CS021HBA1) were flagged by the correlator and 
their data were not used. All 244 SBs were correctly 
processed and stored by the correlator. 

LBA-high (30 — 77 MHz): a second observation was per- 
formed with the LBA system on 14th and 15th of April, 
2011, using a 30 MHz high-pass filter. The visibility 
sampling time was 1 s. The LBA-INNER configuration 
was used. At the end of the data reduction procedure 36 



SBs out of 244 (15%, 7.2 MHz of bandwidth) were not 
usable due to computing-cluster or correlator failure. 
LB A- low (15 — 30 MHz): a third observation was per- 
formed with the LBA system on 16th of July, 2011, 
using a 10 MHz high-pass filter. The visibility sampling 
time was 1 s. Three SBs out of 77 (4%) were corrupted 
during the data processing. We did a visual inspection 
of the 74 residual SBs and only 41 (55%, 8.2 MHz of 
bandwidth) contained usable data, the others were un- 
usable due to high RFI levels. One antenna (CS302) was 
flagged at correlation time. An LBA-OUTER configu- 
ration was used to keep the FoV comparable to that of 
the LBA-high observation. 

International stations were not used in these observations, 
therefore the longest baseline available was about 80 km 
(for the observation at 15 — 30 MHz) and 40 km (for the 
others), while the shortest was 90 m. A plot of the full 
li-y-coverage is shown in Fig. (U while the dirty beams used 
for deconvolution are shown in Fig. [2j 

3.1. The pipeline 

Most of LOFAR data processing is done by pipeline soft- 
ware. The LOFAR processing system takes the data from 
the dipoles, forms the beams and correlates the output to 
ultimately generate images of the radio-sky. As a detailed 
description of the whole process is beyond the scope of this 
paper, here we will illustrate o nly the most import ant steps. 
Interested readers can refer to iHeald et"all (|2010l ) . 

After correlation, data are recorded on storage nodes in 
the current LOFAR offline processing cluster. The first data 
processing step is to flag radio-frequency interference (RFI) 
and optionally compress the data in time and frequency. 
Automated flagging is required for the LOFAR data vol- 
ume, which is performed by the AOFlagger (Offringa e t all 
I2OIOI : lOffringa et aD l2012f ). This software estimates the 
astronomical signal by carrying out a surface fit in the 
time-frequency plane and flags outliers using combinato- 
rial thresholding and morphological detection. Compression 
and further flagging is performed by the New Default Pre- 
processing Pipeline, or NDPPP. 
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Table 1: Details of the observations 



Obs. ID Antenna Frequency Date Observation Sampling FWHM^ Maximum Number 





type 


range 
[MHz] 




length 

N 




time 

M 


[deg] 


resolution 

[arcsec^] 


of stations 


L24923 
L25455 
L29694 


HBA-DUAL^ 
LBA-INNER 
LBA-OUTER 


115 - 162 
30 - 77 
15-30 


2/3-Apr-2011 
14/15-Apr-2011 
16-Jul-2011 


28810 {- 8 
28810 (~ 8 
28805 i~ 8 


h) 
h) 
h) 


1 

2 
2 


- 5 

- 10 

- 10 


19 X 14 
37 X 30 
85 X 44 


45 (7)2 

24 (7)2 

25 (8)2 



Notes. ^ FWHM of the primary beam when points at the zenith, its shape changes during the observation time and is not circular. 
2 "DUAL" means that the two sub-stations of the core stations are treated separately (see text for details). This is why the number 
of stations in the HBA observation is higher with respect to the LBA observations. ^ Enclosed in brackets the number of remote 
stations. 



In cases where the contributions of other bright sources 
in the sky are not negligible, a subtraction of these 
sources directly from the visibilities is required. The tech- 
nique we used is called demixing and it is described in 
Ivan der Tol et all ([2007). This approach is computationally 
cheap and shows remarkably good results in some circum- 
stances, where the source to demix and subtract is not too 
close to the main target, but fails in more complex sce- 
narios, for instance where a strong source is a few degrees 
away from a weak target. We note that the demixing is 
just one of the possible approximate approaches that can 
be used to remove strong interfering sources. Another pos- 
sibility, although computationally more expensive, is the 
peeling procedure (|Noordamll2004l ) . 

The calibration step is performed with the BlackBoard 
Selfcal (BBS) software, developed explicitly for LOFAR. 
This calibration package is based on the Hamaker- 
Bregman-Sault Me asurement Equation ([Hama ker et all 
119961 : ISmirnovll201ll ). which expresses the instrumental re- 
sponse to incoming electromagnetic radiation within the 
framework of a matrix formalism. BBS is thus able to han- 
dle complicated calibration tasks like direction-dependent 
effects and full polarization calibration as well as correcting 
for the time/position dependence of the element beam and 
the synthesized beam. 

The imaging step is routinely performed using the 
AWimager software (Tasse et al. in prep.). The AWim ager 
uses the A-projection algorithm (iBhatnagar et al.]|200 8) to 
image wide fields of view, where data must be corrected 
for direction dependent effects varying in time and fre- 
quency (mainly beam and ionospheric effect s). The software 
does not yet support multi-scale cleaning (|Cornwe III 120081 : 
iRau fc C ornwelll l2QllD . so for our purposes we decided to 
use the CASA[j imager which is adequate for the central 
portion of the image that we were interested in. 

3.2. Data reduction 

Although some steps were common, the data reduction pro- 
cedure followed different schemes for the three observations. 
A preliminary common step is the use of the automated 
AOFlagger on the full resolution (in time and frequency) 
datasets. With a further visual inspection of the data, we 
did not recognise any visible RFI effects left in the raw 
data. After this step the procedures were different for each 
dataset, and therefore we explain them in detail: 



^ http://casa.nrao.edu. 



HBA (115 - 162 MHz): first all baselines related to the 
two substations within the same station were flagged 
(~ 1% of data), this was necessary due to a possi- 
ble cross-talk effect that was found in the intra-station 
baselines. Then, we applied the demixing procedure to 
the dataset, subtracting in this way the two strongest 
sources in the sky, Cassiopeia A (~ 107° from Virgo A) 
and Cygnus A (~ 98° from Virgo A). This procedure 
was necessary only in the second half of the observation, 
where the two aforementioned sources were above the 
horizon. After that we compressed the dataset to one 
channel (excluding the first and last two channels) per 
SB and 20 s of samphng time. This reduced the data 
volume to the level of ~ 400 MB per SB, where a cycle 
of self-calibration lasts ~ 1 h. The model for the self- 
calibration was ex tracted from VLA data at 325 MHz 
(|Owen et al.ll200Qf ), which had a resolution high enough 
for our case (^ 7''). For each SB we rescaled the to- 
tal fiux of the model according to the source global 
spectral index value (see Sect. 13. 3p . Several cycles of 
self-calibration (phase and amplitude) performed with 
BBS and using the imaging algorithm in CASA, were 
necessary to converge to the final image. The imaging 
step at these frequencies was performed using a stan- 
dard CLEAN for the bright central region, followed by 
the use of a multi-scale cleaning. 

LBA-high (30 — 77 MHz): after demixing^ that was per- 
formed as described for the HBA dataset, the data were 
averaged to one channel (excluding the first and last 
two channels) per SB and to 10 s of sampling time. 
The model for self-calibratio n was extracted fr om a 
VLA observation at 74 MHz ( Kassim et allll993f ) with 
the total flux rescaled to the appropriate frequency 
(see Sect. 13. 3p . We did several cycles of self-calibration 
(phase and amplitude) with BBS and using the imag- 
ing algorithm in CASA. The central region of Virgo A 
was CLEANed using standard pixel-by-pixel cleaning 
while for the extended emission we used a multi-scale 
approach. 

LB A- low (15 — 30 MHz): the attempt to use the demixing 
procedure failed since the data taken in the last 3 hours 
of observation were severely affected by ionospheric dis- 
turbances. Since this was the part of the observation 
also corrupted by Cassiopeia A and Cygnus A signals, 
we decided not to use it and simply average the rest of 
data to 5 s and one channel (excluding the first and last 
two channels) before the calibration procedure. Finally, 
several cycles of self-calibration (phase and amplitude), 
using BBS for the calibration and CASA for the imag- 
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Requency [MHz] - LBA Requency [MHz] - HBA 



Fig. 3: Black line: percentage of flagged data. In the frequency range 15 — 30 MHz, the last 37% of the observation was 
manually flagged. Red line: normalized bandpass (for HBA it is only available for a slightly shifted frequency range). 
The SBs removed because of corrupted data or a computer failure are coloured in grey. Completely flagged stations are 
not taken into account to compute the percentage of flagged data. 



ing, were performed. The model for self-calibration was 
again extracted from a VLA observation at 74 MHz with 
the total flux rescaled to the appropriate frequency (see 
Sect. 13. 3p . The imaging step was done in the same way 
as for the LBA-high dataset. At the end of the cali- 
bration procedure, a visual inspection of the images re- 
vealed that for 33 SBs (out of 74) we were unable to 
correctly calibrate the data due to the RFI level. The 
majority of these SBs are indeed concentrated in the 
frequency range 15 — 20 MHz and where the RFI pres- 
ence was critically high. We did not use those SBs for 
the following analyses. 

In Fig. [3] we plot the amount of flagged data for each 
SB, together with the shape of the bandpass functions. The 
amounts of flagged data reflect only partially the amount of 
RFI. Firstly because the RFI flagging is performed at full 
time-frequency resolution and during the subsequent data 
averaging flags are ignored if at least one datum is valid 
in the averaged block. Therefore, it was not possible to 
track those flags due to RFI which are narrow-frequency or 
shorter than the average time. Secondly because new flags 
are applied to remove outliers produced in the calibration 
phase. In the high frequency regime, the percentage of un- 
usable data is more or less constant, at ~ 5%, apart from 
a small increment at 118 MHz. Almost all of these are due 
to the manual-flagging of the first two hours of data from 
RS208 and RS307 and the last two hours from RS208. At 
lower frequencies the RFI is stronger and the peaks in Fig. [3] 
are related to it. Below 30 MHz all SBs had a flagged data 
percentage above 37% because, as explained, we removed 
the last 3 h of observation. In the LBA-high frequency range 
(30 — 77 MHz) the amount of flagged data is rising towards 
the band edges, where it also presents some systematic os- 
cillations. These behaviours are due to the lower sensitiv- 
ity of the instrument at these frequencies, which produces 
some outliers during the calibration procedure and princi- 
pally during the demixing process. These outliers are due 
to a poor signal to noise in the calibration step and were 
flagged by an automated procedure through NDPPP after 
every selfcal cycle. This increases the amount of flagged 
data, but these flags are not RFI-related. The oscillating 
pattern is introduced by flagging outliers after the demix- 
ing procedure, which may suggest that the demixing is less 



effective at those frequencies where the strong {demixed) 
sources are in particular configurations with respect to the 
beam side-lobe pattern. 

3.3. Absolute flux density 

The flux density of Virgo A integrated over all the ex- 
tended emission was rescaled to its expected value to 
compensate for the absence of an absolute flux cal- 
ibrator, while the relative fluxes of different compo- 
nents in the radio morphology is correctly recovered 
by self-calibration. To do that, we collected the total 
flux measurements available in the literatur e in the fre- 
quency range from 10 to 1400 MHz ( Braud e et al.l 119691 
Bridle fc Purtonlll968l:[Ro^er et al.|[T969l: IViner fc EricksonI 



19751; iKellermann et al.l Il969t IWright fc OtruDcekl ll990D. 
Each data-point was corrected to match the iRoger et al.l 
(1973^ (RB C) fl u x sca le with correcti on fa c tors from 
Laing fc PeacocS (|l980r ) and [Scaife fc HeaO (|2012D . A 
model of the form 



log^ = log(Ao) + Ailog 



150 MHz 

V 



(1) 



where v is the observing frequency and S the observed flux, 
was used to flt this data set (see Fig. HJ. 

The model was applied in linear frequency space, i.e. 



N 

= 1 



(2) 



in order to retain Gaussian noise characteristics. 
Parameters were fltted using a Maximum Likelihood 
(ML) approach through a Markov Chain Monte Carlo 
implementation (|Scaife fc Healdl 120121 ). 

We tested polynomial fits up to the fourth order and 
found that a first order polynomial function (Aq = 1226=bl7 
and A\ = —0.79 ± 0.008) is the best-fit model (as alread y 
pointed out by a number of authors, e.g. lRoger et al.iri973f ). 
We derived the expected total flux of Virgo A at the fre- 
quency of each observed LOFAR SB and rescaled the model 
used for that SB to match it at the beginning of each cycle 
of self-calibration. 
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Fig. 4: Integrated flux of Virgo A at different frequencies 
obtained from archival data. The line is a linear fit (slope: 
—0.79) obtained as described in the text. The two vertical 
dashed lines indicate the boundaries of the LOFAR observ- 
ing band. 



The primary beam attenuation at the edge of Virgo A is 
< 3% for HBA images and < 1% for LBA images. We did 
not take this effect into account and included the systematic 
error in the error budget. 

The map at 325 MHz was also rescaled to match the 
RBC flux scale, but at higher frequencies our first order 
polynomial model is probably no longer valid. However, at 
frequencies > 3 00 MHz, the RBC scale i s in agreement with 
the KPW scale ([Kellermann et al. II 19691 ), for which we have 
conversion factors from the Baars scale (|Baars et al.lll977l 
Table 7). Therefore, we used those factors to rescale the 
maps at 1.4, 1.6, and 10.55 GHz from the Baars scale to 
the RBC scale. 



4. Virgo A images 

In Fig. [5] we show the image of Virgo A as seen with the 
LOFAR-HBA at an average frequency of 140 MHz. This 
image is an average over the entire 48 MHz of bandwidth 
and the imaging step has been performed with CASA us- 
ing the multi-scale deconvolution algorithm with a Briggs 
weighting (robust = —0.5). In Fig. [6] we show four images 
of Virgo A obtained with the LOFAR-LBA. Each image 
is realised with CASA usi n^ a multi-scal e multi-frequency 
deconvolution algorithm (C ornwe 1I I2OO8I : iRau fc Cornell 
I2OIID on a subset of 60 SBs (12 MHz of bandwidth) with 
uniform weighting. Finally in Fig. Owe present a very low 
frequency (25 MHz) image of Virgo A. This image was ob- 
tained in CASA using a multi-scale multi-frequency decon- 
volution algorithm with uniform weighting on all usable 
SBs of the LBA-low dataset (20 - 30 MHz). The rms error 
in the images is set by deconvolution errors which limit our 
dynamic range to ~ 5000 (for the HBA map). Our abil- 
ity to recover flux not in the model was confirmed by the 
detection of several sources which were not included in it. 
We recovered > 50 sources in LBA wide-field and > 300 
sources in the HBA wide-field (de Gasperin et al. in prep.). 
The morphological structure of Virgo A instead appear sim- 



ilar from 20 MHz up to the GHz regime, a part from those 
differences caused by spatial changes in the spectral index. 

Virgo A has a ~ 5 kpc-wide inner cocoon, where a one- 
sided jet is visible. The jet, detected also in the optical and 
X-ray bands, points towards the north-west. The counter- 
jet, although not visible due to the effect of relativist ic de- 
boosting, is probably responsible for the emission in the east 
part of the inner region. The inner region is surrounded by 
two much fainter, much larger "bubbles" (~ 40 kpc wide, 
see Fig. [8]) that are overlapping in the central region be- 
cause of projection effects. The inner and the outer haloes 
are connected by two large "flows" (see Fig.[8|). The first is 
oriented almost exactly east- west and the second slightly to 
the north of west, aligned with the inner jet. The eastward- 
flow proceeds straight, forming a well-defined cylinder, and 
ends in a pair of bright lobes, whose edges are brighter 
than their central part. The westward- flow, on the other 
hand, quickly changes its direction projected into the sky 
plane and twists as soon as it leaves the inner region. The 
flow then proceeds towards the south and is compose d of 
a number of thinner structures that, following Owen et al.l 
(|2000t ) , we call "filaments" . Both flows originate in the in- 
ner region and reach the border of the outer haloes. Once 
the halo edge is reached both flows disperse, although only 
the west-flow seems to fill the entire halo with its plasma- 
filaments. The presence of these flows which connect the 
inner halo to the outer edges indicate that the diffuse emis- 
sion is not a simple relic of a previous outburst, but fresh 
energetic particles still flow from the central cocoon. Plasma 
ages derived from spectral fits along the flows, confirm this 
picture (see Sect. 15. 2p . 

The 140 MHz image shown in Fig. [5l although less re- 
so lved than the 7" resolution 325 MHz image presented 
m Uwen et all (|2000D . confirms some characteristics of this 
source. First, the outer halo has a sharp edge and all 
the radio-emitting plasma seems to be confined within 
its boundaries. Second, although less visible than in the 
325 MHz image, part of these edges are limb brightened, 
reinforcing the previous statement. 

Interestingly, we can confirm that this picture is valid 
down to 25 MHz (see Fig. [6] and [7|). If relic emission of past 
AGN activities extending beyond the sharp edges had been 
present, it would have been detectable thanks to its steep 
spectra. However, even at 25 MHz all visible emission is 
confined within the same boundaries that we see at higher 
frequencies. This fact also supports the picture that all of 
the emitting plasma is well confined by the strong pressure 
of the ICM, as we will discuss in more detail in Sect. 15.21 

5. Spectral analysis of the extended halo 

5.1. Spectral index map 

In this section we excluded the very low frequency part 
(^ 45 MHz) of the LBA observations to retain an angu- 
lar resolution of 50''. We produced an image from each 
SB and we convolved them to a resolution of 50''. Then, 
we averaged all LBA and HBA images separately, obtain- 
ing one image in the HBA frequency range and one im- 
age in the LBA frequency range. Finally, the low- frequency 
spectral index map shown in Fig. [9] was obtained by ex- 
tracting a pixel-by-pixel linear regression using three im- 
ages: one extracted from the LOFAR-HBA observation 
(115 - 162 MHz), one extracted from the LOFAR-LBA ob- 
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J2Q00 Right Ascension 

Fig. 5: LOFAR-HBA image of Virgo A at 140 MHz. The rms noise level is a = 20 mJy/beam, the flux peak is 101 Jy/beam 
and the beam size is 21" x 15'' (ellipse in the bottom- left corner). Some deconvolution errors are visible as small holes 
slightly above and below the bright core. 



servation (45 — 77 MHz) and a third one from the VLA at 
325 MHz. 

The central cocoon of the source has a spectral in- 
de:x|3 ranging from —0.55 to —0.6, consistent with what has 
been observed by other authors at higher frequencies up to 
the optical band (Biretta et al. 1991|). The spectrum is a 
straight power-law down to 30 MHz (see Fig. [10]). No ev- 

^ Spectral index definition: oc u^. 



idence of a turnover due to self-absorption is visible down 
to these frequencies. From the total integrated spectrum 
shown in Fig. |4] a possible sign of a turnover in the source 
integrated flux is visible at frequencies < 20 MHz, so out- 
side our frequency coverage. Features north of the bright 
core are likely affected by deconvolution errors and we do 
not consider them as real. 



8 



F. de Gasperin et al.: M87 at metre wavelengths 



(Jy/beam) 



(Jy/beam) 





J2000 Right Ascension 



(c) 59 MHz - RMS: 0.07 Jy/beam - Beam: 45" x 36'' 




J2000 Riglnt Ascension 



(d) 71 MHz - RMS: 0.05 Jy/beam - Beam: 37'' x 30'' 



Fig. 6: LOFAR-LBA images of Virgo A at frequencies ranging from 36 to 71 MHz. Each image is a result of a multi-scale 
multi-frequency cleaning on a subset of 60 SBs. The beam shape is visible in the bottom-left corner of each image. 
Positive contour levels are represented by black lines at (1, 2, 3, 4, 5, 6, 8, 10, 12, 14, 16, 18, 20, 25, 30, 25, 40, 45, 50, 
75, 100) x5<j. Negative contour levels at —la are represented by a red dashed line. 



In the southern lobe the spectral index flattens by 
10% — 20% where the bright flow twists and bends with 
respect to the surrounding areas. The spectral index in the 
east lobe is comparable to what is observed in the southern 
lobe in the radio-brightest zones. The two prominent fila- 
ments to the north-east on the other hand, do not present 



any peculiar spectral index structure, although this can be 
related to the low resolution of our spectral index map. 
The faint extension towards the north-east is the steep- 
est part of the halo, reaching in our map a spectral index 
of —1.8. Interestingly this feature is co- located with what 
iForman et al. (2007) identify as an external cavity in the 
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Fig. 7: Image of Virgo A at 25 MHz. The map noise level 
is (J = 0.6 Jy/beam and the beam size is 85'^ x 44" (grey 
ellipse in the bottom- left corner). Positive contour levels 
are represented by black lines at (1, 2, 3, 4, 5, 6, 8, 10, 12, 
14, 16, 18, 20, 25, 30, 25, 40, 45, 50, 75, 100) x5a. Negative 
contour levels at —la are represented by a red dashed line. 

X-ray halo. The rest of the northern lobe has the lowest sig- 
nal to noise ratio of the map, therefore it could be affected 
by spurious features. 

The northern halo is related to the counter-jet, and 
therefore farther away from the observer. Differences 
between north a nd s outh lobes were also found by 
iRottmann et al.l (|1996[ ) at 10 GHz. They observed a 
higher degree of polarizatio n in the southern lobe (due to 
theLaing-Garrington effect, iLaingifToSSl : iGarrington et al.l 
1 198 81) and a total flux from the southern lobe that is 20% 
higher than from the northern one. 

In general, there is no noteworthy relation between the 
spectral index and surface brightness maps, although a 
steepening of the spectral index is present (at the north- 
east and south-east edges of the map), where a reduction 
in the flow-related activity is present and a flattening is vis- 
ible in some of the flow-active locations (the east lobe and 
the initial part of the west flow). Although in the lowest 
signal-to-noise zone of the map, we report a flattening of 
the spectrum by ~ 20% compared to the rest of the halo 
in the north-west part of it. This feature seems not to be 
related with any structure in the brightness maps. 

5.2. Spectral index fits 

We will now make a detailed analysis of the source radio 
spectrum, in specific regions of interest. For this analysis 
we decided to retain all of the frequencies down to 30 MHz, 
which limits our angular resolution to 75''. The LOFAR 
maps have been averaged in blocks of 10 (bandwidth of 
2 MHz), resulting in 24 maps in the LB A and HBA fre- 



30' 



12°16' 




12^31 "'12'' 06' 00' 30"'54'' 48' 42' 36' 30' 24' 18' 

JZOOO Right Ascension 



Fig. 8: Same as Fig. [5] with most prominent features de- 
scribed in the text labelled. Contour lines are at (1, 2.5, 5, 
7.5, 10, 15, 20, 25, 30, 35, 45, 250, 1000) x3cr. 



quency range each. We also used three archival VLA maps 
at 325, 1400 and 16 00 MHz and a sing;le-d ish Effelsberg 
map at 10.55 GHz (iRottmann et al.l Il996"l ). All of these 
maps were convolved to a resolution of 75 '' and spectra l 
index fits were performed using Synage++ (Murgia 2001"). 
At low- frequencies (< 325 MHz) the spectra are almost 
straight and their slopes are the same reported in Fig. [9l 
although in the following analysis the resolution is lower 
compared to that used for the figure. To assess the reha- 
bility of a spectral index study with images produced by 
different interferometers, the same model of Virgo A was 
simulated in the used LOFAR and VLA configurations and 
frequencies. We imaged those data using the same weight- 
ing scheme (uniform), cell-size and iterations and convolved 
the CLEANed maps at the same resolution to compare the 
outcomes of the different datasets. Virgo A is never resolved 
out, but the different uv-coYeidiges and, to a lesser extent, 
the missing short baselines at higher frequencies, create ar- 
tifacts. Although the ratio between maps produced with 
different instruments shows errors up to 10% in a single 
pixel in the zones where the signal-to-noise ratio is low, we 
note that such errors are not in the form of an overall bias 
but of patches of higher/lower flux (the average error of the 
pixel fluxes across the whole source is +0.3%). However, in 
our analysis we have always used the flux integrated over a 
certain solid angle, therefore we extracted the error for all 
the zones described in Sect. I5.2.2[ finding in every case an 
integrated flux discrepancy below 1%. 

In the following analysis of the spectral data, we tested 
three different models: 

The JP model (|Jaffe fc PerolallT973l ): models spectral age- 
ing as due to synchrotron and inverse Compton losses, 
with the pitch angles of the synchrotron emitting elec- 
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Fig. 9: Left figure: Low-frequency spec- 
tral index map obtained from LOFAR- 
LBA (45 - 71 MHz only) and LOFAR- 
HBA (115 - 162 MHz) observations, to- 
gether with VLA map at 325 MHz. All 
maps were convolved to a resolution of 
50'' (see circle in the lower left cor- 
ner) and a pixel-by-pixel linear regres- 
sion was extracted. Pixels where the er- 
ror were above 3cr are blanked. Contour 
lines are from the 325 MHz map. Top 
figure: spectral index la error map. 
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trons continuously isotropized on a time-scale shorter 
than the rad iative time-scale. 

The KP model (lKardashevlll962l : lPacholczvklll97Q[ ) : as in 
the JP model, but now the pitch angle of the elec- 
trons remains in its initial orientation with respect to 
the magnetic fie ld. 

The CI model ( Pa cholczvklll97Q' ^ : in the "continuous in- 
jection" model, an uninterrupted supply of fresh parti- 
cles is injected by the central source. These particles age 
following the JP model. This model is applicable only 
if the injected particles cannot escape from the selected 
region. Therefore, we used it only on the integrated flux 
from the whole halo. 



Compared to the KP model, the JP model is more realistic 
from a physical point of view, as an anisotropic pitch angle 
distribution will become more isotropic due to scattering, 
magnetic field lines wandering in a turbulent medium, and 
changes in the magnetic field s trength between different re- 
gions (e.g.. ICarilli et al.lll99 1^). Furthermore, since inverse 
Compton losses due to scattering by CMB photons can 
isotropise the electron population, a true KP mo del can be 
visible only for strong magnetic fields (> 30 /iC JSlee et al.l 
I2001D . where inverse Compton losses are negligible. In all 
of these models the injected particles are assumed to have 
a power-law energy spectrum A^(7) cx 7^^^"^ (where 7 is the 
particles' Lorentz factor), which results in a power-law ra- 
diation spectrum with spectral index ainj = {Si^i + 1) /2. 



Finally, the magnetic field strength is assumed constant for 
the entire radiating period. 



Since the most energetic particles radiate their energy 
more efficiently, they are the first to be depleted. Therefore, 
the source radio spectrum evolves and displays a break to a 
steeper slope at a break frequency i^b, which relates to the 
time elapsed from the injection and to the magnetic field 
as B-H-'^ (|Jaffe fc Perolalll973[ ). Major differences 

between the models are visible at frequencies higher than 
i^b, while at lower frequencies all models are expected to 
have a spectral index equal to ainj. There are three free 
parameters in these models: the first is the spectral slope 
ainj of the synchrotron emission generated by an injected 
electron population, the second is the break frequency i^b 
and the third is the overall normalization. 



A modification to the CI model that was also t ested is 
the CIOFF model (Komissarov fc Gubanovl Il994f ). which 
allows the source to switch off after a certain time, encoded 
in an extra free parameter r^off = ^reiic /^source, where ^source 
is the time since the beginning of the outburst and treiic is 
the time since the source switched off. In this case the spec- 
trum would show a first break frequency that depends on 
^source, which Separates two power-laws hke in the standard 
CI model, and an exponential cut-off at higher frequencies 
which depends on treiic- 
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Fig. 10: Fit to the integrated flux spectrum of the central 
region. The dashed line is a linear fit with a slope of a = 
-0.6 ±0.02. 

5.2.1. Central cocoon and macro-regions 

First, a spectral fit was made to the central region (Fig.fTO]) 
defined by C in FiglT3l As the data appear to be described 
by a straight line down to 30 MHz, we fitted a simple power- 
law, obtaining a slope of a = —0.6 ± 0.02. In this and in 
the subsequent fits, the errors on the integrated fluxes are 
computed from the RMS on the individual maps (multi- 
plied by the square root of the number of beams covering 
the area) plus a 10% error due to systematic^ combined in 
quadrature. 

Then, we defined three macro-regions to obtain the av- 
erage spectra of the halo and the flows. 

— The first region (Fig. Illap was obtained by cutting the 
36 MHz map at the Sa level and removing the central 
cocoon. We chose this frequency to maximise the flux 
from the halo (which is higher at lower frequencies). 

— The second region (Fig. Illbp was obtained by remov- 
ing from the previous map all the area with a surface 
brightness above 30 Jy/beam, i.e. the parts of the halo 
dominated by the flows. 

— The final region (Fig. Illcp was obtained by retain- 
ing only the flows-dominated part of the halo (surface 
brightness > 40 Jy/beam), and removing the central 
region. 

A spectral fit using the CI and the CIOFF models has been 
performed on each of these zones and the results are shown 
in Fig. [TTIand in Table [2l Firstly, we performed a standard 
linear regression, from which it can be seen that the spectra 
are curved at high frequencies and, to a lesser extent, also 
at low frequencies. Then, we fitted the data using the CI 
model and fixing ainj = —0.6, equal to the core spectral 
index. In this case the model is not able to reproduce the 
data (see Fig. [TT]) . 

^ The systematic errors considered are: a 1% error due to the 
different 'iz'u- coverages (Sec. 15. 2p . 3% due to the uncorrected 
beam shape (Sec. 13. 3p . and 4% for the uncertainty in the ab- 
solute flux rescaling fSec. 13.3")) . A conservative total systematic 
errors budget of 10% has been adopted. However, this value 
overestimates the real error in the brightest parts of the source, 
therefore providing particularly low Xred values. 
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(a) Halo (without central cocoon) 





CI: i/j,=1.41a=-0,gO 

CI [a fixed]: 1^=0,04 Q=-0.eo 




CIOFF [a fixed]: 0,09 <i=-0.60 Tj^tr =0,04 

POWERLftW: q=-1.15 
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(c) Flows (without central cocoon) 



Fig. 11: Spectral fits of the three macro-regions identified 
in Sect. 15.2.11 In each panel the mask used to select the 
region is shown in black in the bottom- left corner, the white 
part was excluded from the analysis. A simple power-law 
fit (black line) is shown together with two fits of the CI 
and CIOFF models. Red line: fit obtained fixing ainj to 
—0.6. Blue line: with ain} allowed to vary. A CI model with 
ainj = —0.6 is in general not able to fit the data. Finally, 
green lines show the fit of a CIOFF model with ainj fixed 
to —0.6. u\) is in GHz. 
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Table 2: Global spectral fits 



Region CI Model CI Model (ainj = -0.6) CIOFF Model (ainj = -0.6) — Power-law — 

Xred [GHz] -ttinj X^ed [GHz] Xved [GHz] ?7off Xved -^inj 

0.060 0.60t°°^ 

2t'of, 0.0412:^1 1.997 i.nt'of, 

^t'ofs 0.04ig:gl 2.156 1.15ig:g^ 

it'o'ol OMt'o'ol 1-905 1.091^:^2 



Central cocoon 
















Halo (no core) 


0.198 


1 q + 0.3 

^■"^-0.7 




1.471 


U.UO_Q Q]^ 


0.282 


0. 


Halo (no flows) 


0.371 


1 4+0-3 

-■-•^-0.8 


o.9oig:g^e 


2.607 


o.o4ig:gl 


0.439 


0. 


Flows (no core) 


0.121 






0.970 




0.182 


0. 



We decided then to relax some constraints and we re- 
peated the fit using a CI model with all three parameters 
(^b, <^inj and the normalization) free to vary. In this case 
we found a between 1.1 and 1.4 GHz and an ainj be- 
tween —0.83 and —0.9. This result suggests that the ex- 
tended halo is quite young: < 50 Myr, assuming an average 
magnetic field strength of 10 fiG (see Sect. 15. 3p and using 
equation [3l We also observe a steepening of ainj moving 
from the central region (ainj = — 0-6^0 02)5 flows 
(ainj = -0.831^:^7) and the halo (ainj = -O.QOt^;^^), while 
taking the errors into account only a marginal steepening is 
detected moving from the flows to the halo. Although the 
core and the flows are surrounded by the halo, the projec- 
tion effects should not alter these results, in fact the halo is 
on average 20 times and 4 times fainter in flux density 
than the core and the flow regions respectively. In Fig. [12] 
we plot the ratio between the power-law fit to the central 
cocoon data and the CI- model fit to the flow zones, and 
the ratio between the latter and the Cl-model fit to the 
halo without flows. If the emission in these zones were re- 
lated to the same outburst of relativistic particles, simple 
synchrotron ageing would have left the low-frequency part 
of the spectrum untouched at a = —0.6, producing a con- 
stant ratio until the break frequency. Therefore, some other 
mechanism must have steepened the spectra at the lowest 
frequencies. We list here some possible explanations. 

(1) Adiabatic expansion of the relativistic plasma will 
shift the spectra towards lower frequencies and lower in- 
tensities, therefore can also affect the lo w frequency end o f 
the spectrum. A model was dev eloped bv lKardashevI (|1962[ ) 
and revisited by iMurgia et all ([1999). They propose a con- 
tinuous injection of particles which subsequently expand 
adiabatically (CIE model). Such a model produces a low- 
frequency slope of the spectrum that is dominated by the 
plasmas at different ages after adiabatic expansion which 
results in a steepening of the spectra compatible with that 
seen in the lobes of V ir^o A. A similar model was devel- 
oped in lBlundell et al.l ([l999) for double radio sources. The 
authors proposed that adiabatic expansion of plasma may 
happen as soon as the plasma leak out from the hot-spot 
regions. Particles leaving the hot-spots have spent different 
amounts of time in these regions with high magnetic fields, 
experiencing different ageing. Consequently, the final parti- 
cle spectral distribution will be again a sum of spectra with 
many break frequencies and adiabatically expanded. This 
idea was claimed by Blundell & Rawlings (2000) to explain 
the "i njection index discrepancy" discussed in Carilli et al. 
(|l99lf ). In this last paper, the authors observe that the low 
frequency spectral index measured in the lobe of Cygnus A 
should reflect either the low-frequency spectral index of the 
hot-spot {a = —0.5, marginal adiabatic losses) or the high- 
frequency spectral index of the hot-spot {a = —1, strong 
adiabatic losses). They observe instead a low- frequency 
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Fig. 12: The solid line shows the ratio between the power- 
law fit to the cocoon zone (Fig. [10]) and the Cl-model fit to 
the flow zones (Fig. Illcp . The dashed line shows the ratio 
between the Cl-model fits to the flow zones (Fig. Illcp and 
to the halo (Fig. Illbp . Simple synchrotron ageing would 
have left untouched the low-frequency part of the spectrum, 
producing in this plot a horizontal line until the point where 
the break frequency occurs. We notice instead a steepening 
in the spectrum, going from the central region to the flows 
and, to a lesser extent, from the flows to the halo. 



spectral index of —0.7. These numbers are not particularly 
different from what we observe in Virgo A. Although these 
sources are remarkably different, the underlying physical ef- 
fect which dominates at the low-frequency end of the spec- 
trum might be similar. 

(2) During their lifetime, relativistic electrons crossed a 
wide range of magnetic field strengths. As a result, their 
final spectrum is the sum of many spectra with many dif- 
ferent break frequencies. Particles that have spent much of 
their life within strong magnetic fields will have a very low 
break frequency, which can modify the low frequency slope 
of the radio spectrum. 

(3) A third scenario is that the radio spectrum is intrin - 
sically curved even in the core (iBlundell Rawlingsll2QQQl) . 
In the source core the magnetic field strength is few mG 
while it is o:^ 10 /iG in the halo (see Sect. 15. 3p , implying 
that the radio emission from the core is powered by elec- 
trons with energy ~ 20 times smaller than those emitting 

1/2 

in the lobes, according to j oc {u/B) ' ,7 being the typical 
Lorentz factor of electrons emitting at the frequency in a 
magnetic field of strength B. Thus, it is possible that the 
low-frequency synchrotron spectrum of the extended halo 
reflects a curvature of the spectrum of the emitting particles 
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in the core at higher energies; this is expected in the con - 
text of particle acceleration models (Amato fc Aronsir2QQ6f ). 
Fig. [To] does not provide compelling evidence in favour of 
this scenario, as the spectrum does not appear to steepen 
significantly at higher frequencies, yet future observations 
with higher resolution will test this possibility. 

(4) Another scenario is that we are observing the relic 
emission of a source that had been active for an extended 
period and recently stopped injecting plasma in the halo. In 
this case the flatter spectral index (ainj — —0.6) tail of the 
spectrum lies at lower frequencies, where we cannot detect 
it, and the ain} — —0.85 slope is the steepened part of the 
spectrum. Following a simple CI prescription, we should 
expect the high-frequency part to have a slope of —1.1. 
However, to fit the data, a simple CI model is not enough 
and an exponential cut-off must occur at a frequency of ~ 
5 GHz, implying a recent switching off of the fresh particles 
flow. This model (CIOFF) retains an initial injection index 
of ainj = —0.6, but fits the data much better than the 
standard Cl-model with the same initial slope (see Fig. [TTI 
and Table [2]). In the CIOFF case the halo is older (thaio — 
150 Myr, assuming an average magnetic field of 10 /iC and 
using equation [3]) than in the CI scenario and the source 
must have switched off only a few Myr ago. This "switching 
off" should be interpreted as the most recent detaching of a 
bubble from the source central region, while in the central 
cocoon a new bubble is forming. In this case, the steeper 
low-frequency slope is the consequence of the integration 
over a wide area of the source, which again implies the 
addition of many spectra with different break frequencies 
as a consequence of ages differences. 

5.2.2. Zone by zone spectral analysis 

To check the validity of these hypotheses through an in- 
depth spectral analysis, we selected ten relevant zones (de- 
fined in Fig. [13]) considered to be representative of the dif- 
ferent parts of the source: central cocoon, fiows/filaments 
(zones Wl, W2, W3 and W4 for the west flow; zones El 
and E2 for the east flow), and halo (zones HI, H2 and H3). 
We deliberately avoided those zones where the signal to 
noise ratio is low and data could be dominated by calibra- 
tion or deconvolution errors. For each zone we extracted 
the averaged flux at each frequency and performed a spec- 
tral analysis (see results in Table [3]). Since our data are 
well approximated by a straight line in the LOFAR bands, 
with the curvature determined only by the three higher fre- 
quency values, both JP and KP models are able to fit the 
data comparably well (see Fig. [14]) when leaving free to 
vary all the parameters. To reduce the degrees of freedom, 
we fit the data by fixing ainj = — 0.6 and ainj = —0.85 (see 
Fig. [IS]). 

A fixed ainj = —0.6 fails to fit our data, and when all 
of the parameters are left free to vary, we find z^b ^ 5 GHz. 
These results are in contradiction with the CIOFF scenario 
(model 4 in Sect. I5.2.ip . where we would have expected 
a wide range of z^b and ainj = —0.6. On the other hand, 
we can fit both these and the global spectra if we assume 
an ainj — —0.85. This steeper slope can be due to the in- 
trinsically different and rather steep injection spectrum of 
a previous outburst. However, in that case, the presence 
of uninterrupted flows from the cocoon to the outer halo 
would be harder to explain. 



(jy/beam) 

5 10 15 20 




12^3nO^ 00' 30^50^ 40' 30= 20^ 

J2000 Right Ascension 



Fig. 13: Image of Virgo A at 31 MHz convolved to a resolu- 
tion of 75". The ten different zones we analysed have been 
outHned in the map. 

Alternatively, we argue that a global steepening of the 
spectrum may occur at a very early stage, when the bubble 
detaches from the central cocoon. This can be related to an 
adiabatic expansion of the bubble, which happens as soon 
as it l eaves the high-pr essure central region urazov et al. 
200li ICarilh et al.l [l991) . In this case, a mix of plasmas 
at different ages would quickly expand and the high fre- 
quency parts of t heir spectra, w hich is indeed curved and 
steeper than —0.6 (jCotton et al.l l2Q09). would be shifted to- 
wards lower frequencies. The sum of these spectra produces 
a resulting spectrum that has a low- freq uency end steeper 
than the initial ones (iMurgia et al.lll999[ ). Furthermore, an 
abrupt lowering in the spectra normalization is also ex- 
pected, this can account for the strong brightness contrast 
between the cocoon and the halo. After that no further 
strong expansions occurred, otherwise a gradient in the sur- 
face brightness across the halo would be visible. 

Regardless of the mechanism responsible for the ob- 
served steepening, in what follows we assume ainj = —0.85 
at the point where the plasma bubbles leave the cocoon. 
This provides a thaio — 40 Myr (from equation [3l based on 
the break frequency in the CI model fit to the entire halo 
and assuming an average magnetic field strength of 10 jaG). 

5.3. Magnetic fields and synchrotron ageing 

Equipartition magnetic field strengths were computed in 
the same regions shown in Fig. [131 We made the standard 
assumption that each zone of the radio source contains rel- 
ativistic particles and uniformly distributed magnetic fields 
under energy equipartition conditions (e.g. Miley 1980). 
The radio luminosities for the computation were extracted 
from the lowest frequency map (31 MHz). To perform this 
analysis we made the following assumptions: 
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Table 3: Spectral fits to representative regions 



Region JP Model KP Model — Power-law ■ 







Xred 


i^h [GHz] 


Q^inj 


Xred 


i^h [GHz] 


Q^inj 




' Wl 


0.090 


15.6t^i 




0.083 


O q + 0.3 

0.0_2 5 




West ^ 


W2 


0.028 






0.029 


^•^-1.2 


0.8llH^ 


Flow 


W3 


0.060 






0.047 


4.8i?-J 


0.87in^ 




< W4 


0.034 






0.033 




0.891H? 


East 


f El 


0.022 


7.6t;-j 


0.83lH^ 


0.010 


o 4+0.2 


PI O-I+0.02 


Flow 


[ E2 


0.034 






0.042 


4 Q+0-9 
^•^-1.3 






f 


0.056 


8.3+1-? 


0.82lH^ 


0.068 


4 1+0-2 


0.82lH^ 


Halo < 


H2 


0.034 


-'--'-•^-2.6 


0.82lH^ 


0.035 


6.0+;-^ 


0.821H? 




. H3 


0.045 


Q 2+0-9 


o.83in^ 


0.044 


4 6+°-^ 


0.82in^ 



2 _ 
Xred ^inj 

1.551 1.04t°;°^ 



3.030 



3+0.02 



3.120 1.19lgg^ 

3.259 1.23tH^ 

4.370 1.191°;°} 

2.943 1.231°-°^ 

3.088 1.15t°;°^ 

2.229 1.09t°;°2 

2.997 1.141°;°^ 



Notes. The assumed 10% error in the flux densities overestimates the real random, normally distributed error providing artificially 
low x?ed values. 



The particle energy U is equally divided between elec- 
trons and protons, setting Upr/Ue\ = k = 1. However, 
in Table |4] we also list the results for the assumption of 
/c = 0, where all of the particle energy is in the relativis- 
tic electrons (and positrons), in line with the idea that 
the jets could be main ly composed of these constituents 
([Revnolds et aljri996[ ). We show also the results for the 
assumption of A: = 10, as un upper limit from what pre- 
dicted by electromagnetic acceleration models involving 
protons and electrons (Bell 1978). 

The equipartition magnetic field is usually computed 
assuming that the relativistic particle energies are 
confined between a minimum emin and a maximum 
Cmax value, corresponding to an observable frequency 
ran^e, typically a ssumed to be 10 MHz - 100 GHz 
(|Pacholczvklll970[ ). However, a fixed frequency range 
corresponds to an energy range that depends on the 
magnetic field, which varies in different parts of the 
source. Therefore, we decided to put limits directly 
on the electron populatio n energies ([Brunetti et al.l 
Il997l : [Beck fc Krausd [20051 ) . This approach, compared 
to the standard one, provides slightly higher B values 
(see A ppendix A, Brunetti et al. 1997). Reynolds et al. 
(|l996l ) and iDunn et al.l ()2QQ6f ) put constraints on the 
maximum value of 7min noting that the synchrotron 
self-Compton flux density generated in the source core, 
which depends on 7min, cannot exceed the observed X- 
ray flux density. They obtained 1 ^ 7min < 100 for 
an electron-positro n jet and 50 < 7min 5, 1 00 for an 



electron-proton jet. iFalcke fc Biermaiml (|l995l) also ar- 
gue for radio loud AGN a 7min — 100. We repeated our 
analysis for two values of 7min* 10 (with k = 0) and 100 
(with k = 1 and with k = 10), corresponding to emin 
of 5 and 50 MeV respectively. Above 7min — 1000 we 
would expect a turnover in the low frequency part of 
our spectra that we do not detect. The emax value does 
not affect the results, and we used an arbitrarily high 
value of 5 GeV (7max ^ 10000). 

For each zone we assumed a cylindrical configuration 
and we repeated the computation for two depth D = 
20 kpc and D = 40 kpc. In the rest of the paper the 
flow zones are assumed to have a depth of 20 kpc and 
the halo zones are assumed to have a depth of 40 kpc. 
The depth of the core is assumed to be 5 kpc. 
The low frequency spectrum slopes have been assumed 
to be equal to —0.85, as observed. 



The equipartition magnetic fields, the minimum pres- 
sures and the corresponding synchrotron ageing times are 
listed in Table [H In the first part, we list the equiparti- 
tion analysis results assuming that all of the energy re- 
sides in the electrons and positrons only {k = 0) and that 
7min = 10. In the middle and third parts of the table, we 
assume a 7min values of 100 and we relax the electrons to 
protons energy ratio to A: = 1 and A: = 10 respectively. In 
the rest of the paper we will refer to the values in the second 
part {k = 1 and 7min = 100) only, while we notice that the 
other two combinations of parameters provide ages so small 
that the plasma bubbles should move at a velocity higher 
than the sound speed (c^ 900 km s~^ in the outskirts). 
The source age is obtained by (see e.g. iMurgia et al.l 

^0.5 

ts = 1590 (3) 

where the synchrotron age is in Myr, the magnetic field 
strength in /iG and the break frequency in GHz, while 
Bic = 3.25(1 + z)^ /iG is the inverse Compton equivalent 
magnetic fiel d strength with energy density equal to that 
of the CMB (jSlee et al.ll2"QQll ). The break frequencies were 
obtained from the fit of the JP model (see Table [3|). We as- 
sumed a constant and uniform magnetic field strength and 
neglected any influence on the spectra from e.g. expansion 
or local re-energization of electrons. 

The equipartition analysis provides us with reasonable 
values for the lifetimes of the bubbles. For example, follow- 
ing the various zones sampled in the west flow (Wl, W2, 
W3, and W4), a bubble escapes fro m the source cocoon 
after 7 Myr (|Churazov et al.ll200lL estimated - 10 Myr 
with simulations) and reaches the outer edge of the lobe 
after 12 — 15 Myr (zone W4). This confirms the picture 
of a flow continuously replenished with fresh particles. In 
the east flow the centre of the lobe (zone E2) is reached 
after 10 — 13 Myr. This lifetime is about a fourth of what 
is derived by dynamical models. A difference between these 
age estimations and the global halo age {c^ 40 Myr) is ex- 
pected, as the latter takes into account many regions which 
have a lower break frequency than those in the flow zones. 

In Fig. [16] we plot the theoretical temporal evolution 
of the radio spectrum in the halo zone HI using the stan- 
dard JP model. We simulated how the spectrum of a zone 
with these characteristics can evolve and how old it must 
be in order to go undetected in our maps. We find that 
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JP: 1-^=15.58 a=-0.83 



JP: 1/^=9.06 a=-0.82 





-- JP; 1/^=9.79 a=-0. 88 











(a) Zone Wl 



(b) Zone W2 



(c) Zone W3 



JP: i'i,=7.57 Q=-0.83 



JP; 1/^=9.62 a=-0.91 



(d) Zone W4 



(e) Zone El 



(f ) Zone E2 



JP; 1-^=8. 34 a=-0.82 



X 



JP; 1/^=11.44 a=-0. 82 



X 



JP; i'i,=9.24 a=-0.83 



(g) Zone HI 



(h) Zone H2 



(i) Zone H3 



Fig. 14: Fit of the JP (blue), KP (red) models to the zones related to the halo. The zones are defined in Fig. [131 The 
black line is a simple linear regression fit to emphasize the curvature in the spectrum. u\) is in GHz. 



in this illustrative example we would be able to detect 
emission as old as ~ 400 Myr. This number is rather op- 
timistic and it would decrease if adiabatic expansion had 
played an important role or if we had overestimated the 
current age of the zone. Finally, given the confinement of 
the source discussed in Sect. 13. 3[ it is possible that parti- 
cles from older AGN events were mixed with those from 
recent events. This can also play a role in steepening the 
low-frequency end of the lobes' spectrum. Such a scenario 
was observed in simulation by Morsony et al. (2010), where 
the authors simulated AGN driven jets in a dynamic, cos- 
mologically evolved galaxy cluster. They found that largest 
scale reached by AGN jets is only proportional to the AGN 

1 /3 

power {R ex Pj ) and does not depend on the activity 
time. In this case the estimated halo age should be inter- 
preted as a lower limit. 



6. Discussion 

Radio morphological evidence (Sect. H]) and spectral anal- 
ysis (Sect. [5]), show that the Virgo A halo is an active part 
of the source and not a relic of past activities (as already 
pointed out bv*Owen et al.l l20QQl ). Its radio emission is con- 
fined inside sharp boundaries. 

The low-frequency spectral index map is fairly uniform, 
apart from a flattening in the central cocoon and in the 
northern lobe, and a steepening in the regions where the 
flow activity is fading. Thanks to the LOFAR data, together 
with high frequency observations up to 10 GHz, we were 
able to extract wide-band radio spectra of the source halo 
to an unprecedented detail. A continuous injection model 
applied to the whole halo shows cut-off frequencies at ^ 
1.3 GHz, which provides an estimation of the halo age of 
40 Myr. 

We performed a detailed spectral analysis of nine dif- 
ferent zones in the halo. Leaving all the parameters free 
to vary, we obtain a good fit only assuming ainj — —0.85 
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Fig. 15: Fit of the JP (dotted cyan and magenta hnes) and KP (dashed blue and red lines) models for the zones related 
to the halo. The zones are defined in Fig. [131 The slope of the injected electron population was fixed to ainj = —0.6 (cyan 
and blue lines) or ainj = —0.85 (magenta and red lines). The black line is a simple linear regression fit to emphasize the 
curvature in the spectrum, z^b is in GHz. 



(see Fig. [Hj). In order to reduce the number of free pa- 
rameters we tried to set the slope of the energy distri- 
bution of the injected electron population to 5in} = —2.2 
(<^inj = —0.6), as we observe in the source central region, 
and which is also in perfect a greement with the broad - 
band spectrum of the jet of M87 (jPerlman fc Wilson|[2QQ5l[ ) . 
However, in this case the JP and KP models fail to repro- 
duce the observed spectra, while fixing ainj — —0.85 pro- 
duces reasonable results with both models (see Fig. [T5]) . 
Observationally, this reflects the presence of a steep low- 
frequency end in the spectra of all regions in the halo, even 
those just outside the central cocoon. We speculate that this 
steep spectrum is the consequence of the strong adiabatic 
expansion of a mix of plasmas at different ages that takes 
place as soon as the plasma bubbles leave the dense cen- 
tral area, shifting their already steepened high-frequency 
part of their spectrum down to the MHz region. Once we 
assume ainj — —0.85, we are indeed able to follow the age- 



ing of the plasma bubbles along their path. For example, 
in the west flow we observe a shift of the break frequency 
(15 GHz 9 GHz 8 GHz), which reflects an ageing of 
the electrons (7 Myr 10 Myr 12 Myr), as we move far- 
ther from the source centre along the flow. It is worth noting 
that between region W2 and W3 the results with error-bars 
are consistent with no-ageing, we can speculate that those 
two regions are, in projection, at the a similar distance from 
Wl but on different fronts of the flow, which is then creating 
a mushroom-like structure close to the edge of the halo (as 
in the east flow) instead of flowing along it. This picture can 
also explain the spectral index map in the southern lobe, 
where the flatter values, which seem to be connected with 
the active flow and initially follow it, end directly at the 
edge of the lobe that is where the flow should separate and 
create the top of the mushroom-shaped lobe. This picture is 
very much in line with the source three dimensional recon- 
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Table 4: Equipartition analysis 
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Notes. D is the depth of the region assuming a cylindrical configuration. Beq and pmin are the magnetic field and the pressure 
from the equipartition analysis, t is the estimated zone age. Errors on t are derived from errors on z^b in Table [3] 



struction proposed in IChurazov et al] (|200lf ), but requires 
more detailed analyses to be confirmed. 

The minimum-pressure analysis in nine different zones 
in the source provides equipartition magnetic field values 
and relative pr essures. The se values are compatible with 
those found by I Owen et all ((20001), although some of the 
initial assumptions were different: a different low-energy 
cut-off, different spectral slope and a different model for 
the synchrotron emissivity. Nevertheless, the net balance 
of these modifications reduces only slightly the values of 
the magnetic field strength we estimated for the Virgo A 
halo. It is worth pointing out that all the obtained val- 
ues are a few times above 3.25 /iG, below which the Inverse 
Compton scattering of CMB photons is the dominant effect 
compared to synchrotron losses. 

We compared our minimum p ressure values with thos e 
derived from XMM observations (' Matsushita et all [20021 ). 
listed in the last column of Table HI The thermal gas pres- 
sure is always more than an order of magnitude higher than 
the minimum pressure obtained from equipartition analy- 
sis. This could imply that the regions where plasma elec- 
trons and magnetic fields are co-located are much smaller 
than the assumed zone volumes (a "filling factor" < 1), or 
that much of the energy is contained in the thermal gas, 
or that the source is not at the equipartition. Yet another 
possibility is that much of the energy is stored in the rela- 
tivistic protons (see e.g. lPunn fc Fabianll2004l : fBirzan et al.l 
I2008D . In this case, to keep the equipartition principle valid, 
we would have to assume a > 1000, although this would 
generate much higher magnetic field strengths and therefore 
unrealistic short plasma ages. It is well known that inside 
the inner cocoon radio jets create cavities in the thermal 
gas, whereas on the larger halo scales there is evidence for 
the uplift of cold t hermal gas (see e.g. iForman et al.ll2QQ7l : 
iMillion et al.ll2QlQl ) along the radio-flows. This has been ex- 



plained by IChurazov et al.l (|2QQlf ) as gas uplifted by buoy- 
ant bubbles from the cold cluster core. In such a scenario 
the amount of pressure provided by thermal gas cannot be 
neglected and accounts for a not-negligible fraction of the 
total pressure. Furthe rmore, using X-ray and optical data, 
IChurazov et al.l (|2QQ8[ ) also found that the combined contri- 
bution of cosmic rays, magnetic fields and micro- turbulence 
to the total pressure in the core of M87 is ~ 10% of the gas 
thermal pressure. This is in line with our findings. 

6.1. Plasma age and dynamical time 

Hydro-dynamic simulations of buoyant bubbles suggest 
that the halo o utskirts can be reached in 42 — 67 Myr 
(Churazo v et al. lf2QQl, the proper number depending on the 
orientation and assumed distance). These estimations are 
a factor of ~ 4 greater tha n what was found in our syn- 
chrotron spectral analysis. iGull fc NorthoveJ (|l973t ) first 
noticed that the lifetime of synchrotron-emitting particles 
is short compared with dynamical time-scales. Some rea- 
sons that can e xplain this discrepancy are p resen ted in 
'Ch urazov etlIlj2QQlf ). iBlundell fc Rawlingsl (|2QQQD . and 
Owe n et all (|2QQQ[ ): (i) the bubbles may be filled with a mix 
of weak and strong magnetic fields, the relativist ic electrons 
may survive for a long time in the weak magnetic fields and 
then radiate most of their energies as soon as they diffuse in 
the strong magnetic field regions. This picture would also 
account for the filamentary structure visible in the radio 
images, (ii) I n situ acc eleration of particles may play an im- 
portant role (iParma et al. 1999 : Prieto et al. 2002) . In both 
these cases the age of the radiating particles can hardly be 
determined by radio observations, (iii) New plasma can also 
flow along pre-existing channels and replenish aged plasma. 

Finally, the buoyant rise time may not necessarily be 
the correct time, the radio plasma that forms bubbles 
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Fig. 16: To estimate our capability to detect ancient out- 
bursts we determined the maximum age of an emission we 
would be able to detect given the sensitivity of our maps. 
We took the zone HI (Fig. [T3|) as representative for an old 
and inactive zone in the halo. Using its = 8.3 GHz and its 
estimated age of 18.7 Myr we let its spectrum evolve. In this 
figure we show the fit to the observed data as a black line. 
Coloured lines show the evolution of the expected spectral 
shape as time passes by (from top to bottom, after 250, 300, 
350, 400, 450, 500 and 550 Myr, respectively, after the elec- 
tron injection). Triangles indicates the 3<j detection limit 
for our maps. Blue is for the LB A- low map (Fig. [7]), red for 
the four LBA-high maps (Fig. [6]), while yellow for the HBA 
map (Fig. [5|). 



comes from a jet whi ch likely leaves it with a momentum. 
iBriiggen et al.l (|2002l ) simulated this situation in an envi- 
ronment compatible with that of Virgo and found in this 
case that the plasma can reach the distance of 20 kpc in 
:^ 15 Myr, which is very much similar to what we measure. 



6.2. Energetics 

In I Owen et al.l (|2000( ) the authors derived an estimate of the 
halo age using energetic conservation arguments, following 
lEilek fc Shorl (119891) . from 



dt 



dV 



^rad- 



(4) 



where p the pressure of the halo plasma, V the volume, Uint 
the internal energy, Pj the jet power and Lrad is the radia- 
tive losses from the halo, which is small compared to Pj and 
can be neglected. We can solve the equation assuming that 
the halo is spherically symmetric with a radius R and that it 
expands due to its own internal energy up to R = 35 kpc. 
We also assumed that Pj is time-independent. If the ex- 
pansion is slow, the pressure of the bubble can be approxi- 
mated with the pressure of the surrounding medium, which 



we o btained from XMM t emperature and density pro- 
files (Matsushita et al"ll2002[ ). For a non-relativistic plasma 
(most of the bubble content is thermal, with F = 5/3), we 

obtain a halo age of - 250 (P^/10'^'^) Myr. For a bub- 
ble dominated by relativistic particles and magnetic fields 

(F = 4/3), we obtain - 400 (P^/IO"^"^) Myr. It is im- 
portant to notice that a shorter burst of energetic particles 
would require a much higher Pj (although for a shorter 
duration) than a long continuous injection of particles. In 
this simple model, our measurements of the halo age based 
on equipartition {t :^ 40 Myr), provide an estimate of the 
jet power of Pj 6 x 10^^ erg s~^ for F = 5/3 and of 
Pj 10 X 10^^ erg s-i for F = 4/3. This result is con- 
sistent with the conclusions of lOwen et al.l (^000) where 
the jet po wer Pj is estimated to be ^ few x 10^^ erg s~^ 
and with iDi Matteo et al.l (|2003t ) which predict a Pj 
5 X 10^^ erg s"-*^ for accretion at the Bondi rate. The re- 
sult is ins tead above wha t has been found by other au- 
thors, e.g. iRevnolds etld] (|1996D find Pj 10"^^ erg s-^ 
This energy supply is ^ 100 times higher than the X-ray 
luminosity of the cooling flow region around M87, that is 
- 3 X 10"^^ erg s"^ (Mushotzky & Smith 1980). Therefore, 
even a modest efficiency of energy dissipation into heat, 
is able to exceed the radiative cooling of the gas. Finally, 
we notice that an energetically proton-dominated plasma 
{k 1000) would lower by a factor of ~ 10 the age of 
the halo, boosting the necessary jet power up by the same 
amount which would drive the source at quasar luminosity 
(|Falcke fc Bier mam] 1 19951 ), which is not observed. 

6.3. Ultra high energy cosmic rays 

A debated argument is the possibility for a radio galaxy 
like M87 to accelerate ultra high energy cosmic rays 
(UHECR) of ^ lO^Q eV di rectly in the radio lobes. The 
iPierre Auger Collaboration! (|2007[ ) found a correlation be- 
tween the arrival directions of cosmic rays with energy 
above 6 x 10^^ eV and the positions of AGN within ~ 
75 Mpc, with a small excess of detections in the direction 
of Centa urus A, but no e vent detected in the direction of 
Virgo A. IChisellini et al.l ()2008h correlated the position of 
UHECR events with the directions of Hl-selected galaxies, 
and proposed that the UHECR coming from the direction 
of Centaurus A instead originate from the more distant 
Centaurus cluster, whose galaxies are in fact richer in Hi 
than the galaxies of the Virgo cluster, explaining why there 
is no UHECR event from the Virgo direction. This opens 
the possibility that UHECR are produced by gamma ray 
bursts or newly born magnetars, but does not rule out an 
AGN origin. 

In Centaurus A, the lobe radius 100 kpc wide) 
and estimated magnetic field strength {B 1 /iC, 
iHardcastle et al.l 2009) satisfy the Hillas argument (Hilla^ 
"1984'), which requires that the Larmor radius of the ac- 
celerated particles fits within the source. The maximum 
particle energy Ep = £^20 10^^ eV that can be accelerated 
inside the Virgo A lobes is 
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Z e tl B Z rioo B-e eV , 



where rioo = 100 kpc and B-q = B 10~^ G. Using an 
average magnetic field of 10 /iC and a lobe radius of 20 kpc, 
we estimate the maximum particle energies to be Ep 
Z 2 X 10^^ eV, where Z is the atomic number. However, 



19 



F. de Gasperin et al.: M87 at metre wavelengths 



efficient stochastic acceleration processes due to resonant 
interaction between particles and magnetic field turbulent 
perturbations require lobe contents to be predominantly 
relativistic plasma, and this is probably not the case for 
Virgo A. 

7. Conclusions 

In this paper we presented the first LOFAR observation 
of Virgo A in the 15 — 162 MHz frequency range and an 
analysis of its radio spectra. The major results are: 

— Virgo A halo is an active part of the source and not 
a relic of past a ctivities (as already pointed out by 
lOwen et aLll2QQQf ). 

— Going down to 25 MHz no previously unseen steep- 
spectrum features were detected. The source appears 
instead to be well confined within boundaries that are 
identical at all frequencies. 

— A low- frequency spectral index map of Virgo A shows no 
obvious relation between spectral index and brightness. 
Generally, a steepening of the spectral index is present 
where the morphological evidences of flow activities are 
reduced in the southern lobe and at the position of the 
north-east X-ray cavity. A flattening of the spectral in- 
dex, instead, is visible at the northern lobe towards west 
and in the position of enhanced radio flux density. 

— A spectral analysis of the extended halo was performed. 
With the assumption of ainj = —0.6 as in the core re- 
gion, characterized by the flatter spectral slope, none 
of models tested is able to flt the data. Instead, ainj = 
—0.85 is required, outside the central cocoon, to explain 
the spectra. This steepening in the low-frequency end of 
the spectra can be connected to a strong adiabatic ex- 
pansion of the plasma bubbles that happens as soon as 
they leave the dense central region. 

— An equipartition analysis was conducted and an aver- 
age magnetic field strength of :^ 13 /iG was found in 
the flow regions, while a magnetic field of 10 /iG is 
present in the halo regions. In the inner cocoon the av- 
erage magnetic field reaches 30 /iG. 

— A synchrotron ageing analysis provided a global halo 
age of ^ 40 Myr. The particle age tends to increase 
with distance from the centre in the fiow regions and 
reaches a maximum of 15 Myr where the fiows end. 
This age is about a factor of ~ 4 less than the dynamic 
time of a buoyantly raising bubble. 

— Minimum pressure analysis reveals that, given our as- 
sumptions, the pressure generated by the plasma and 
the magnetic fields is much less than what is required 
to sustain the halo against external pressure. Probably, 
thermal gas is also playing a role in sustaining the halo. 
On the other hand, some of the parameters we assumed, 
or the equipartition hypothesis itself, may not be cor- 
rect. 

— We estimate the jet power Pj to be 6 — 10 x 10^^ erg s~^. 
This energy supply is 10 to 100 times higher than the 
X-ray luminosity of the cooling fiow region around M87 
(~ 3 X 10^^ erg s~^) and is therefore able to exceed the 
radiative cooling of the gas. 

The extended radio-halo of Virgo A is a composite of 
many plasma bubbles of different ages. They are infiated 
in the central cocoon, where clear cavities in the thermal 



gas are visible, by the source's powerful jets. As soon as 
they leave the dense central area, a strong adiabatic ex- 
pansion shifts their radio emission towards lower frequen- 
cies and lower fiuxes. Then they buoyantly rise towards the 
halo outskirts where they disperse. During their motion, 
the travelling bubbles lift up cold X-ray emitting gas from 
the centre of the cluster. 

To better constrain different synchrotron models and 
ultimately draw conclusions about the underlying mecha- 
nisms, further observations are necessary. Steps forward in 
the study of Virgo A require: (i) high angular resolution 
maps of the Virgo A halo at frequencies above 4 GHz, (ii) 
higher angular resolution and (iii) higher dynamic range 
on low frequency (< 200 MHz) maps. Once these require- 
ments are fulfilled, wide-frequency spectral studies can be 
performed on regions of the halo that are small enough to 
disentangle the diffuse emission from that originating in the 
filaments. 

Point (i) will be soon feasible with the VLA and next- 
generation single-dish radio telescopes, like the Sardinia 
Radio Telescopes. In the C-band (4-8 GHz) the VLA 
will require mosaicing and single dish observations to fill in 
the shortest baselines, the results of these observations will 
be presented in forthcoming papers. At higher frequencies, 
single dish radio-telescopes have the required sensitivity to 
detect the extended halo and they can be used at the price 
of lower resolution, currently ~ 40'' at 20 GHz. 

Requirements (ii) and (iii) will be soon available with 
new LOFAR data and presented in forthcoming papers. 
The present antenna configuration (excluding international 
stations) can provide baselines of 80 km, which allow an 
angular resolution of 25'' at 30 MHz. Furthermore, the im- 
portant technical improvements achieved during commis- 
sioning will provide better knowledge of the telescope be- 
haviour while the 8 new stations now available will give a 
much better uv-coYeidige. With these improvements we will 
be able to greatly reduce deconvolution errors and improve 
the dynamic range of future maps. 
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